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Abstract We provide new exact solutions to the Einstein-Maxwell system of 
equations which are physically reasonable. The spacetime is static and spheri- 
cally symmetric with a charged matter distribution. We utilise an equation of 
^ ' state which is quadratic relating the radial pressure to the energy density. Earlier 

models, with linear and quadratic equations of state, are shown to be contained 
in our general class of solutions. The new solutions to the Einstein-Maxwell are 
found in terms of elementary functions. A physical analysis of the matter and 
electromagnetic variables indicates that the model is well behaved and regular. In 
particular there is no singularity in the proper charge density at the stellar centre 
unlike earlier anisotropic models in the presence of the electromagnetic field. 
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The study of charged relativistic objects in general relativity is achieved by solv- 
ing the Einstein-Maxwell system of equations and imposing conditions for physical 
acceptability. This is not easy to achieve because of the nonlinearity of the field 
equations. The exact solutions found have many applications in relativistic astro- 
physics. The models generated have been used in the description of neutron stars 
and black hole formation by Ray et al. [1] and de Felice et al. [2]. Particular models 
have also helped in the establishment of the absolute stability limit for charged 
k> ' spheres by Giuliani et al. 3_ and Bohmer and Harko ^ . Several models of charged 

5_j , relativistic matter have been used to study strange stars by Mak and Harko [5], 

' Bombaci [S], Komathiraj and Maharaj [7] and Thirukkanesh and Maharaj [H]. 
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Charged models have been also used in the description of strange quark matter 
by Discus et al. [9], hybrid protoneutron stars by Nicotra et al. [10], and bare 
quark stars by Usov et al. [llj . A geometric approach is to assume the existence 
of a group of conformal motions on spacetime; exact solutions have been found 
by Mak and Harko |12j for strange quark matter and Esculpi and Aloma [H] for 
anisotropic relativistic charged matter by assuming the existence of a conformal 
killing vector in static spherically symmetric spacetimes. 

Models with an equation of state are desirable in the description of realis- 
tic astrophysical matter. However most explicit solutions of the Einstein-Maxwell 
system that have been found do not satisfy this property. There have been some 
attempts made recently to find exact analytic solutions of the Einstein- Maxwell 
system with a linear equation of state. These include the treatments of Ivanov [TJ , 
Sharma and Maharaj [15] , and Thirukkanesh and Maharaj . Particular solutions 
with a quadratic equation of state, relating the radial pressure to the energy, where 
found by Feroze and Siddiqui [16] . This is an important advance since the com- 
plexity of the model is greatly increased because of the nonlinearity of the radial 
pressure in terms of the energy density. However the investigations mentioned 
above all suffer from the undesirable property of possessing a singularity in the 
property charge density at the centre of sphere. An essential requirement for a well 
behaved electromagnetic field is regularity of the proper charge density through- 
out the matter distribution, particularly at the stellar centre. The importance of 
this feature has been highlighted in the analysis of Varela et al. whose treat- 
ment offers a general approach of dealing with anisotropic charged matter with 
linear or nonlinear equations of state. It is desirable to eliminate the singularity 
in the charge density for a detailed and complete analysis of physical properties 
of charged compact objects. 

Our results may be helpful in the study of compact stars and gravitational col- 
lapse relating to neutron stars and black holes. In this regard we refer to particular 
papers some of which have static spherical geometry and others are dynamical. 
Novikov [18] showed in the case of spherical geometry that collapse of electrically 
charged matter may replaced by expansion and infinite densities are avoided. A 
general treatment of collapsing charged matter was completed by Bekenstein [19] 
who showed that nonzero pressure plays a significant role. The analysis of Ray- 
chaudhuri [20] for charged dust distributions showed that conditions for collapse 
and oscillation depend on the ratio of matter density to charge density. If this ratio 
is large, corresponding to weakly charged dust spheres, then shell crossings cannot 
be avoided in gravitational collapse as proved by Ori [21] . Krasinski and Bolejko 
[22| showed that there exist initial conditions for a charged dust sphere with finite 
radius so that a full cycle of pulsation can be completed by the outer layer with 
no internal singularity. A full and comprehensive analysis of charged, dissipative 
collapse is provided by Di Prisco et al ^23^ for the free-streaming and diffusion 
approximations. A related and detailed analysis in the gravitational collapse of a 
charged medium was performed by Kouretsis and Tsagas [24 where the role of 
Raychaudhuri equation is highlighted. Exact solutions with an equation of state, 
such the quadratic case considered in this paper, are helpful in such studies. 

The objective of this paper is to find new exact solutions of the Einstein- 
Maxwell field equations with a charged anisotropic matter distribution and a 
quadratic equation of state. We indicate that particular models found in the past 
with an equation of state are part of our general analytical framework. Previous 
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solutions with a linear or quadratic equation of state are regained in our treatment. 
We ensure that the charge density is regular at the centre of the compact body 
and the physical criteria are satisfied. In Sect. [51 we give the Einstein-Maxwell 
field equations for a static spherically symmetric line element as an equivalent 
system of differential equations utilizing a transformation due to Durgapal and 
Bannerji ^5\. In Sect.|3l we present new exact solutions to the Einstein-Maxwell 
system with a quadratic equation of state. The solution is regular at the centre 
of the compact object. This analysis extends the treatment of Thirukkanesh and 
Maharaj [8], and Feroze and Siddiqui [16]. Known solutions with an equation of 
state are presented in Sect. 21 as particular cases of our new results. In Sect. [5l a 
physical analysis of the new solutions is performed; the matter variables and the 
electromagnetic quantities are plotted. We summarise the results obtained in this 
paper. 



2 Field equations 

In standard coordinates the line element for a static spherically symmetric space- 
time, modelling the interior of the relativistic object, has the form 

We take the energy momentum tensor to be of the form 

T,, = diag(-p - pr - pt + pt + ifi"), (2) 

where the quantity pt is the tangential pressure, pr is the radial pressure, p is the 
density, and E is the electric field intensity. Then the Einstein-Maxwell equations 
can be written in the form 

^Hl-e-'')]' = P+lE^ (3) 
-^(l-e--) + ^e--=p.-ii^^ (4) 

-2xf II , 12 , v' , I A'\ 1 2 /t-s 

e \v +v H V \ =pt-|- -E , (5) 

V r r / 2 

- = y2^-\^"E)'^ (6) 

where primes represent differentiation with respect to r. The quantity a represents 
the proper charge density. We are utilising units where the coupling constant 
= 1 and the speed of light c = 1. The mass within a radius r of the sphere is 

M{r) = ^£ u^piu)cLj. (7) 

We now introduce a new independent variable x and define new functions y and 
Z so that 

x = Cr\ Z{x) = e-^^^'"\ ^VW = e2-('"), (8) 



4 



S. D. Maharaj, P. Mafa Takisa 



where A and C are constants. We assume an equation of state of the general form 
Pr = Pr{p) for the matter distribution. We take the quadratic form 



Pr 



7P +ap- (3, 



(9) 



relating the radial pressure pr to the energy density p. In the above a, and 7 
are constants. Then the Einstein-Maxwell equations governing the gravitational 
behaviour of a charged anisotropic sphere, with a quadratic equation of state, are 
represented by 



P_ 
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(10) 

(11) 
(12) 



(13) 

(14) 
(15) 



where dots denote differentiation with respect to the variable x. Equations pOI) - 
psp are similar to the field equations of Thirukkanesh and Maharaj [8_ ; however 
in this case the equation of state is quadratic. The quantity A = pt — pr is called 
the measure of anisotropy and vanishes for isotropic pressures. The nonlinear sys- 
tem as given in (|10p - (|15p consists of six independent equations in six variables 
involving the matter and electromagnetic quantities p,Pr,Pt, A, E,a and the two 
gravitational potentials y and Z. The nonlinearity of the Einstein-Maxwell system 
(|10p - (fT5| ) has been increased, when compared with many earlier treatments, be- 
come of the appearance of the quadratic term in (|lip : when 7 = then there is 
a linear equation of state. In addition equation p4p now contains terms with E'^ 
increasing the complexity of system since 7 7^ in general. 



3 New solutions 



To integrate the Einstein-Maxwell system we make the particular choices 

1 + ax 
E'^ _ k{3 + ax) + sd^x^ 



G 



(1 + ' 



(16) 
(17) 
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The gravitational potential Z is well behaved and finite at the origin. The electric 
field intensity E is continuous, regular at the origin and approaches a constant 
value for increasing values of x. The constants a, &, k and s are real. The general 
analytic functional forms for Z and E regain particular cases studied in the past 
with an equation of state. 

On substituting (|16p and (|17p into (|14p we obtain the first order equation 



ij _ (l + a)(a-6) a(a-h) /3(l + aa;) 

y ~ 4[1 + (a - h)x\ 2(1 + aa;)[l + fex] ~ 4C[1 + 6a;] 

(1 + Q)[fc(3 + ax) + sa^x^\ C7[(3 + ax)(2a - ab - k) - sa^x^f 



(1 + ax)\l + hx\ 



+ ■ 



16(l + aa:)(l + 6a;) 



,(18) 



for the metric function y. In spite of complexity of equation (|18|) it can be solved 
in general. On integrating p8|) we get 



2/ = 15(1 + aa;)"'[l + 6x]" exp {F(x)\ , 



(19) 



where D is the constant of integration. The function F(x) is given explicitly by 



F(x) = 7[2(a- 6) - k] 



2(26 - a)(l + ax) + (6 - a) 



8(6-a)2(l + aa;)2 



-Cs7 



(a - 6)^(aa; + 2) - a(2a + s)(l + i 



4(a-6)(l + aa;)2 
(a - 6)(4fc + s) + (2fc(6 - 3a) + 36s)(l + ax) 



32(a - 6)2(1 + ax)2 



(20) 



The constants m and n have the form 



(l + a)(s + 2fc) a ,i2 

m = --^^ — ^ '- H 1-7 2(a- 6) - fcr 

8(6- a) 2 ^ ' ' 



+ 



1 

+ - 



(6-a)3 (6-a)2 ' 4 



, [(a - 6)[2s(a - 6) + a + 6] + 3a6(fc - 26) - 6^fc 

+26^(2a - 1)] 

n=il±i^l[2(a-6)-fc] + i-l + '^)^-2-('^-^) ' " ^) 



86 



+7[2(a-6)-fc]' 



4(6- a) 



a)3 (6-a)2 



+ 



4C62 

Csa2(l + a) 
862(6- a) 



Cs7 



1662(6 -a): 



[a^(s + 46) + (fc + 26)(6a^6^ - 2a^6)] 
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Then we can generate an exact model for the Einstein-MaxweU system (|10p -(ll5 p 
in the form 



2A 1 + ax 



P_ 

C 

Pr 
Pt 

A 
C 



a 
~C 



: A^D^ (1 + axf"" [1 + bxf' exp[2F(a;)], 

(2a - 26 - fc)(3 + ax) - sa'^x^ 
2(1 + aa;)2 ' 

: 7p^ + ap - /3, 
■Pr + A, 

4:x{l + bx) m{m — l)a^ , 2mnab 



(21) 

(22) 

(23) 

(24) 
(25) 



1 + ax 



{l + axy 



^ 2maF{x) ^ 



{l + ax){l + bx) 1 + ax 



b^n{n - 1) 2nbF{x) 
(1 + bxy ^ 1 + bx 

bn 



+ F{x) + F{x) 



+ 



2{a-b)x 4(1 + 5a;) 
■(l + aa;)2 + 



(1 + ax) 



1 + aa; ^ 1 + 6a; 
1 



+ F{2 



C(2(a - 6) - fc)(3 + ax) - Csa'^x'^ 



2(1 + 1 



2C7(a- 6) + fc(3 + ax) + , 



2(l + aa;)2 
- -^(^ + ax) ] 



(2(a - 6) - fc)(3 + ax) - sa^x'^ 



2(l + ax)2 



E'^ _ k{3 + ax) + sa^ 

IT ~ 



C[l + bx] l^Vk{a^x^ + 3aa; + 6) + 2^50x^3 + ax{2 + ax) 
x{3 + ax){l + ax)^ 



(26) 
(27) 

(28) 



where F(x) is given by (|20p . We observe that the exact solution (|2ip - (f28| ) of the 
Einstein-Maxwell system has been written solely in terms of elementary functions. 
For this solution the mass function is given by 



M{x) 



[(12a(a - 6) - 6ak)x + s(15 + lOaa; - 2a'^x'^)]x^^'^ 



8C3/2 

5 s arctan(Y^ax) 



3a(l + ax) 



7,3/2 



(29) 



The gravitational potentials, matter variables and electromagnetic variables are 
well behaved and regular in the stellar interior. However in general there is a sin- 
gularity in the charge density at the centre which is evident in (|28p . This singularity 
is avoidable when fc = 0, so that we have 



_ 4Csa^a;[l + 6a:](2 + ax] 
"C " (l + aa;)5 



(30) 



At the centre of the star a; = and the charge density vanishes. 
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4 Known solutions 



We have found a general class of exact solutions to the Einstein-Maxwell system 
with a quadratic equation of state. It is interesting to observe that for particular 
parameter values we can regain uncharged anisotropic and isotropic models (fc = 0, 
s = 0) from our general solution (f2T |) - (f28|) . We regain the following particular cases 
of physical interest. 



4.1 Feroze and Siddiqui model 

This is a special case of our general solution with the quadratic equation of state 
Pr = 7p^ + ap — 13. If we set s = 0, C = 1 and A^D^ = B, then we regain the line 
element 



ds 



B{1 + ar^ni + br^r exp[2F(r)]d<2 + 



dr 



+r'^{d9'^ + sin^ 



(31) 



where 



F{x) = -^+^[2{a-b)-kf 



Ab 

a (1 + a)k 
2 ~ 4(6 - a) 

+ a 



2{2b - a){l + ar'^) + {b - a) 



?(6-a)2(l + ar2)2 



+ 7[2(a-6) -fc]' 



(6-a)3 {b-ay 



[2(a - 6) - fc] + ii±i^fc-Hfi(^ + - ^) 



-7[2(a- 6) - kY 



4(6 -a) 



(6-a)3 (6-a)2 4 



462 



The line element (|3ip was found by Feroze and Siddiqui [T^ which was the first 
model with quadratic equation of state. Some minor misprints in [TB] have been 
corrected in our result. This solution may be used to model a compact body. 



4.2 Thirukkanesh and Maharaj model 

If we set 7 = then we have the linear equation of state pr = ap — /3. Also setting 
C = 1, s = and b = a — b, we get the line element 



ds^ = A'D^il + ar'^f"'[l + (a - 6)r^]^" exp 



2(a- 6) 



dr 



+ - 



1 + ar^ 



1 + (a - b)r' 



-dr^ + r^{d6^ + sin^ 



(32) 
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where 

2ab- {l + a)k 



m 



Ah 

n = — [2a^(fc(l + a) - 2ab) - a&(5fc(l + a) - 26(1 + 5q)) 

+6^(3fc(l + q) - 26(1 + 3q) + 2/3)] . 

The metric p2p was found by Thirukkanesh and Maharaj [8]. This solution may 
be used to model realistic charged compact spheres and strange stars with quark 
matter in the presence of electromagnetic field. 



4.3 Sharma and Maharaj model 

If we set 7 = 0, /3 = ap, then we regain the linear equation of state pr = a{p — p) 
where p is the density at the surface. By setting fc = 0s = 0,6 = a — 6, C = l and 
we find the following form of the line element 

, 2 D^i 1 2N2mr, , / 7n 2i2n ( — a/3r^ \ 2 

ds = —B[l + ar ) [1 + [a — b)r \ exp ^ dt 



2{a-b) 



1 + ar^ 



H ' = dr +r(d9^ + sin edd)^), (33) 

l + (a-6)r2 



where 



a 

m = — , 
2' 



5aba - 2a^a - ib^a + a6 - 6^ + 6/3 

n = = . 

4(a- 6)2 

The line element ([33)) represents an uncharged anisotropic sphere and was found 
by Sharma and Maharaj [15] . It may be used to describe strange stars with a linear 
equation of state with quark matter. 



4.4 Lobo model 



If we set 7 = 0, /3 = 0, then we obtain the linear equation of state pr = ap. On 
setting k 
element 



setting fc = 0, s = 0, 6 = a — 6, a = 26, C = 1 and A^D^ = B we regain the line 



where 



ds^ = _(1 + 26r2)2'"(l + br^f^dt' + ( ) dr^ 

\^ 1 + 6r2 J 

+r'^{de'^ +sin ed(l?). (34) 



a 
2' 



1 - a 
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The metric p4(l was first found by Lobo [26] which represents uncharged anisotropic 
matter. This solution serves as a stehar interior with a < — ^ and may be matched 
to the Schwarzchild exterior for dark energy stars. 



4C(1 -x)' 



4.5 Isotropic models 

We observe that Z\ ^ in general and the model remains anisotropic. However, 
we can show for particular parameter values that Z\ = in the general solution 
([21]) -([281). If we set b = {a - 1), k = 0, s = 0, a = 0, then we obtain 

a 

m = — 
2 

3(l + a)C][/3- (1 + 3q)C]. (35) 

Two different cases arise from (I35p by setting A = 0. Firstly, we observe that when 
P = and a = — 1 then A = 0. The equation of state becomes Pr{= Pt) = —p- 
With the line element 

ds^ = - + j dt^ + ^1 + j dr^ + r^(d0^+sin2 0d0^) (36) 

where we have set A = D = \ and C = We mention that the metric (|36p is 
the isotropic uncharged de Sitter model. Secondly, we observe that when /3 = 
and a = — \ then Z\ = 0. The equation of state becomes Pr(= pt) = —\p with the 
line element 

ds^ = -A^dt^ + ^1 + ^ j dr^ +r^(d6>^ +sin^6ld<?i^) (37) 

where D = \ and C = The metric (13711 is the isotropic uncharged Einstein 
model. 



5 Physical Analysis 

We show that the exact solutions of the Einstein-Maxwell system found in section 
3 are well behaved by generating graphical plots of matter and electromagnetic 
variables. We used the software package Mathematica [27 and we make the par- 
ticular choices C = 1, a = 2.5, 6 = 2, 7 = 0.01, a = 0.33, s = 0.017 and k = f3 = 0. 
We generated the plots for the energy density (Fig. 1), radial pressure (Fig. 2), 
electric field intensity (Fig. 3), charge density (Fig. 4), mass (Fig. 5), speed of 
sound (Fig. 6), tangential pressure (Fig. 7) and the measure of anisotropy (Fig. 
8). The energy density p is a finite and monotonically decreasing function. The 
radial pressure pr is similarly well behaved and continuous. The electric field in- 
tensity E is initially small and approaches a maximum value as the boundary is 
approached. The proper charge density a is nonsingular at the origin, increases 
and then decreases after reaching a maximum value. The mass function is a strictly 
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0.5 1.0 1.5 2.0 25 3.0 

Fig. 2 Radial pressure 



increasing function which is continuous and finite. The speed of sound is less than 
the speed of light and causality is maintained throughout the stellar interior. The 
radial pressure is decreasing and does reach a finite value of the radial coordinate. 
The tangential pressure is also a decreasing function. The measure of anisotropy 
is a decreasing function as the boundary is approached and remains finite in the 
interior. Thus all the matter variables, electromagnetic variables and the gravita- 
tional potentials are nonsingular and regular in the region containing the stellar 
centre. In particular the proper charge density a is finite at the centre unlike earlier 
treatments. 

It is desirable to study comprehensively the stability of our new models; this 
is a objective for the future reseach. The solutions generated may be matched to 
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Fig. 5 Mass 




Fig. 8 Measure of anisotropy 
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the exterior Reissner-Nordstrom spacetime 



1 r / \ r r'^ j 

+r^{de'^ +sin^ ed(t>'^), (38) 
across the boundary r = 5R. This generates the following conditions 

2M 

and 



-1 



l-^ + C = ^V(CSR^ (39) 



2M ^l + aC^ (40) 

^ W j 1 + bCW ' ^ 

relating the constants a, b, A, C, a, P and 7. There are sufficients number of free 
parameters to ensure the continuity of the metric coefficients across the boundary 
of the star, ft is possible to study the astrophysical significance of the exact so- 
lutions to the Einstein-Maxwell equations found in this paper. This is the object 
of future research. We point out that for suitable parameter values we regain the 
mass M = 1.433Mq of Dey et al [28 -[30_ corresponding to a strange star model 
when there is no electromagnetic field. Therefore the solutions found in this paper 
may be used to generalise earlier results and to model charged relativistic strange 
and quark stars. 

Our aim in this paper was to find new regular exact solutions to the Einstein- 
Maxwell system for spherically symmetric gravitational field with an eq of state. 
In particular we selected a quadratic equation of state relating the energy density 
to the radial pressure. The new models presented in this paper may be used to 
model relativistic compact objects in astrophysics. 
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